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In this work the dynamics of a chain consisting of a set of beads attached to the 
ends of segments of hxed lengths is investigated. The chain fluctuates at constant 
temperature in a viscous medium. For simplicity, all interactions among the beads 
have been switched off and the number of spatial dimensions has been limited to 
two. 

In the limit in which the chain becomes a continuous system, its behavior may 
be described by a path integral, in which the rigid constraints coming from the in- 
finitesimally small segments are imposed by means of a functional delta-function. In 
this way a model of the dynamics of the chain is obtained, which closely resembles a 
two-dimensional nonlinear sigma model. The partition function of this generalized 
nonlinear sigma model is computed explicitly for a ring-shaped chain in the semi- 
classical approximation. The behavior of the chain at both long and short scales of 
time and distances is investigated. The connection between the generalized nonlinear 
sigma model presented here and the Rouse model is discussed. 



I. INTRODUCTION 

Subject of this work is a study of the dynamics of a continuous chain which is subjected 
to thermodynamic fluctuations at constant temperature T. The chain is represented as the 
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limit of a system of beads and links of fixed length, in which the number of beads becomes 
infinite, while the length a of the links goes to zero in such a way that the total length of the 
chain L = Na remains a finite constant. Problems of this kind are encountered for example in 
polymer physics, because with some approximation long flexible polymers may be regarded 
as continuous chains [Ij. It is thus spontaneous to consider an isolated polymer fluctuating 
in some viscous environment as a concrete realization of the system investigated here. A 
few applications in which the dynamics of a chain turns out to be relevant are mentioned in 
Ref . ^ . In particular, in Ref . six regimes of the chain dynamics are distinguished, which 
apply not only to the well known case of a polymer in a solution, but also to other cases, 
like for instance those of an isolated cold chain and of a hot polymer in the vapor phase. 

From polymer physics, one may borrow the standard approach to the dynamics of a 
chain. It consists in considering the fluctuations of the chain as a stochastic process, which 
is usually described with the help of the Langevin equations or, alternatively, of the Fokker- 
J'lanck-Smoluchowski equations This approach leads to the well known models of Rouse 
31] and Zimm 4[] which allow a satisfactory understanding of the main properties of polymers 
in solutions. One major drawback of these coarse grained models is that they do not take 
into account the rigid constraints which are necessary in order to keep constant the length 
of the links connecting the beads. The Rouse and Rouse-Zimm equation consider instead 
beads joined together by springs, where the local spring is infinitely extensible. In this 
way the length of the chain is not fixed and is allowed to become infinite. Moreover, in 
the continuous limit the Rouse equation is nothing but the stochastic equation (Langevin 
equatiori) for the classical Wiener measure, which yields paths without well defined tangent 
vectors 1]. These problems have been tackled by various attempts, see e.g. js, 6|. However, 
the correct use of rigid constraints in (stochastic) dynamics requires some mathematical 
effort 0, S], in contrast to the static cases where rigid constraints can be implemented by 
Dirac delta functions in the partition functions. For instance, the probability distribution 
of ideal closed chains topologically linked together may be represented as a path integral 
of noninteracting particles with the insertion of Dirac delta functions which enforce the 

"n — 



topological constraints on their trajectories 9|, [lOl, lUl, Il2 |. 

Here a strategy similar to that used for static chains will be applied to dynamics. We 
consider the distribution "^disc of the probability that a fluctuating chain passes from a initial 
discrete spatial conformation Fj to a final one Fj in a given interval of time At = tf — ti. 
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The idea behind our approach is based on the fact that the beads of the chain may be 
regarded as a set of Brownian particles with constrained trajectories. The constraints 
arise due to the presence of the — 1 hnks of fixed lengths connecting the beads. As 
a consequence, it is possible to write the probability distribution "^disc in the form of a 
path integral describing the fluctuations of A^ Brownian particles with the insertion of Dirac 
delta functions. The latter are needed in order to impose the necessary conditions on the 
trajectories of the particles. For simplicity, possible interactions among the beads have been 
switched off, including hydrodynamic interactions and only the two-dimensional case has 
been discussed. 

The limit in which the chain becomes continuous is not entirely trivial. It involves the 
vanishing of three crucial quantities, the mass of the beads, their mobility and their size. 
After performing this limit very carefully, we obtain as a final result the probability distri- 
bution \E' of the continuous chain. It is found that \E' consists of a path integral which closely 
resembles the partition function of a two-dimensional nonlinear sigma model 13|. For this 
reason the obtained model will be called the generalized nonlinear sigma model (GNLSM). 
The main difference from the nonlinear sigma model is that the holonomic constraint is 
replaced in the GNLSM by an nonholonomic constraint, which requires that the tangent at 
every point of the trajectory of the chain is an unit vector. The Lagrange multiplier that 
imposes this nonholonomic constraint plays the same role of pressure in the hydrodynamics 
of incompressible fluids, a fact already noted in Ref. 21]. The "incompressibility" is related 
to the fact that it is not possible to "compress" the lengths of the links joining the beads. 
The GNLSM is well suited to study all situations in which it is possible to assume that the 
thermodynamic fluctuations and the conformational changes are small and slow. Systems 
satisfying these requirements are for instance cold chains or chains in a very viscous medium. 

In principle, the GNLSM is exactly solvable after performing two gaussian integrations, 
but the presence of non-trivial boundary conditions complicates the calculation of the proba- 
bility distribution \E'. Even the method of the effective potential, which is useful to investigate 
phase transitions in nonlinear sigma models, cannot be easily applied. As a matter of fact, 
in the nonlinear sigma model the effective potential is computed assuming that the field 
configurations which minimize the action are constants. However, in the GNLSM configura- 
tions of this kind correspond to the situation in which the chain has collapsed to a point and 
are thus unphysical. Despite these difficulties, it is possible to compute the probability dis- 
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tribution \1/ and the associated generating functional of the correlation functions of the bond 
vectors "^[J] using a background field approximation, in which small gaussian fluctuations 
are considered in the background of a dominating classical conformation. The initial and 
final conformations of the chain are picked up by choosing in a suitable way the background 
classical solution and by tuning the boundary conditions of the gaussian fluctuations. The 
semiclassical approximation is valid in the case of low temperatures or of highly viscous 
media, exactly the regimes in which the GNLSM can be applied. The explicit formulas of 
\1/ and "^[J] derived here show that the fluctuations which deform the chain along directions 
which are tangent to the trajectory of the classical background conformation propagate dif- 
ferently from the normal fluctuations. This fact is a direct consequence of the presence 
of rigid constraints and could be relevant for example in the theory of formation of single 
polymer crystals 1^ . 



The material presented in this paper is organized as follows. In Section II the form of 
the Lagrangian of a classical discrete chain in two dimensions has been derived in polar and 
in cartesian coordinates. No restrictions are posed to the motion of the discrete chain. The 
limit to a continuous chain is however performed assuming that one of the ends of the chain 
is fixed. The calculation in polar coordinates shows that one crucial term in the Lagrangian 
disappears in the continuous limit. This fact simplifies the classical equations of motion of 
the chain. The probability distribution \E' of the fluctuating chain is computed in Section III 
using a path integral approach. Some subtleties arising when taking the continuous limit 
in the probability distribution of the discrete chain are discussed. Section IV is dedicated 
to the study of the classical solutions of the GNLSM. It is shown that the only possible 
classical solutions are time independent, apart from rigid translations with constant velocity 
of the whole chain. The computation of the probability distribution \1/ and of the generating 
functional "^[J] for a ring-shaped chain is performed in the semiclassical approximation in 
Section V. The physical interpretation of the results obtained in the previous Sections is 
presented in Section VI. The equilibrium limit of the GNLSM and its connection with the 
Rouse model are studied in Section VII. Finally, in Section VIII the conclusions are drawn 
and possible future developments are discussed. 
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II. THE CLASSICAL DYNAMICS OF A CONTINUOUS CHAIN 

Let us consider a discrete chain of — 1 segments of fixed lengths I2, . . . ,In in the two- 
dimensional plane. Each segment Pi+iPi is completely specified by the positions of its end 
points Pj+i and Pi. In cartesian coordinates (x, y) these positions are given by the radius 
vectors: 

Yii = {x,,yi) i = l,...,N (1) 

The segments are joined together at the points Pj, where 2 < i < N — 1, see Fig. [1], while 
Pi and P/v are the ends of the chain. Moreover, at each point Pj, with i = 1, . . . , A^, a mass 
irii is attached. In the following we restrict ourselves to the case of a free chain. We will see 
below that the addition of interactions is straightforward. 

At this point we compute the kinetic energy of the above system: 

N 

= E V (*^' + y^) (2) 

i=l ^ 

The subscript disc in the left hand side of Eq. ([2]) is to recall that we are considering at the 
moment a discrete chain with — 1 segments. For future purposes, it will be convenient to 
introduce the kinetic energy of point Pf. 

K. = ^ix^ + y^) (3) 

Of course, to Eq. ([2]) one should also add the constraints 

(x, - + (yi - y^.^f = i = 2, . . . , N (4) 

in order to enforce the requirement that the segments have fixed length /j. It is possible to 
eliminate these constraints passing to polar coordinates: 

i i 

^i = Yl h COS v^j 2/i = E h sin V^j i = l,...,N (5) 

j=i k=i 

ifj is the angle formed by segment j with the y— axis, see Fig. [H According to our settings, 
the radial coordinates Ij are constants for j = 2, . . . , A^. The length li, which denotes the 
distance of the end point Xi, yi from the origin, is not fixed, so that li = li{t) is allowed to 
vary with the time t. From Eq. ([5]) the velocity components of the z— th segment may be 



written as follows: 



= — ^ Iji^j sin ipj — liipi sin ipi + li cos ipi 

i-l 

?/i = ^ COS ipj + /ivji COS ifi + /i sin 



z = 2,...,iV 



(6) 



z = 2,...,iV 



(7) 

xi = — sin yji + /i cos v?! (8) 
= /iv^i cos v^i + /i sin y9i (9) 

Separating the contribution coming from the first i — 1 variables as shown in Eqs. (jSHTj), 





FIG. 1: A chain with TV segments. Let us note that the end point Pi is not bound to stay at a fixed distance 
with respect to the origin of the cartesian reference system. 



the kinetic energy Ki of the i— th point can be expressed in terms of the kinetic energy Ki^i 
of the {i — 1)— th segment: 

i-l 



Ki = —^Ki_i + -^if^pf + miY^ klj^i^j cos{(pj - Lpi) + rriiliipili sin(v9i - Lp,i) (10) 

It is possible to solve the above recursion relation and to find a closed expression of Ki in 
polar coordinates. If we do that, at the end the total kinetic energy of the discrete chain 
becomes: 

M - ^ 



Kh,. 



disc 



= -TT {kVi + ^i) + ^iV^i H cos{(fi-j+i - (fi) 

^ i=\ j=\ 

N i-l 

+ h'YYl "mik-j+iLpi-j+i sin(v9i - v^i-j+i) 

i=l 3=1 
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N i-l N i-1 i-j 

+ li-j+i^'ftj+i + rriik.j+ilkVi-j+i'pk cos((^i_j+i - ^k) (11) 

2=1 j=l 1=1 j=l k=2 

where M = X^iLi is the total mass of the chain {2^. 

We wish now to perform the hmit in which the chain of — 1 segments becomes a 



continuous system 



24l |. To this purpose, it is convenient to consider the indices k, . . . 



appearing in Eq. (fTTj) as discrete variables taking values in a one-dimensional lattice with 
points. Quantities /, carrying the index i may be interpreted as discrete functions of i. Their 
variations A/j are given by: A/j = /j+i — /«. Clearly, Ai = 1, i. e. the spacing between 
two neighboring points in the lattice is 1. In order to proceed, we rescale the distances in 
the lattice in such a way that the spacing in the new lattice will be a. To this purpose, we 
perform the transformations i — > Si, fi — > f{si) where the new variable Si has variation 
Asi = Si+i — Si = a. The next step is to compute the kinetic energy of Eq. (fTTj) in the limit 

N — ^00 a — >0 Na = L (12) 

in which the product Na remains finite and is equal to the total length of the chain L. 
Clearly, in the limit f|T2l) the right hand side of Eq. ffTTl) will diverge unless we suppose 
that the masses rrii and the lengths k of the segments are going to zero in a suitable way. 
Reasonable assumptions are: 

k — > l{si) = acr{si) rrii — > m{si) = ap{si) (13) 

where cr(sj) and p(sj) are respectively the distribution of length and of mass along the chain. 
To be consistent with our settings, the distributions a{si) and p{si) must be normalized as 
follows: 

N N 

Y,a{si)As, = L P{si)^s, = M (14) 

i=l i=l 

While it would be interesting to study chains in which the segments are allowed to have 
different lengths and the points have different masses, for simplicity we will suppose from 
now on that the length and mass distributions in the chains are uniform, i. e.: 

M 

a(si) = 1 for i = 2, . . . , AT and p{si) = — foTi = l,...,N (15) 

Ij 

In the discrete case (compare with Eq. (fT3!l ) this implies that all segments of the chain and 
the masses rrii are equal: 

M 

h = a ioT i = 2, . . . , N and rui = —a for i = 1, . . . , N (16) 

Ij 
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At this point we are ready to pass to the continuous hmit. Functions of discrete variables 
will be substituted with functions of continuous variables, while sums will be replaced with 
integrals according to the following rules: 

f{s,)-^f{s) E^^^^ [ (17) 

i=l 







After a few calculations one finds: 

Kdiscit) Kit) (18) 



where 



Kit) = fimm+im 

+ <pi{t)li{t)— I ds dui^{t, s — u) cos{ip{t, s — u) — Lpi{t)) 



L Jo Jo 

+ li{t)— I ds I dui^it, s — u) sm{{pi(t) — {pit, s — u)) 
L Jo Jo 

/ ds du dvip{t, s — u)Lp{t,v) cos{Lp{t, s — u) — (p{t,v)) (19) 
L Jo Jo Jo 

Let us note that the right hand side of Eq. f|T9|) contains four terms, while the original discrete 
version of the kinetic energy in Eq. (ITT]) contained five terms. In fact, the contributions 
proportional to ^f-j+i of Eq. ( fTTI) disappear in the continuous limit. 

Eq. ( fT9l) may be simplified by performing the following change of variables: 



u' = s — u du' = —du (20) 

Using also the formula: 

ds r du'fiu') = r ds(L - s)f(s) (21) 
JO Jo Jo 

which is valid for any integrable function f{s), we obtain: 

M 

+ '^i{t)hit)-r ds{L - s){p(t,s)cos{ip{t,s) - ipiit)) 



L Jo 
M 

+ kit)— ds(L - s)0(t, s) sin(v9i(t) - (/9(t, s)) 
L Jo 

+ — / ds{L — s) / duipit, s)ip{t,u) cos{ipit, s) — (fit,u)) (22) 
L Jo Jo 
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As a further simplification, one could fix the point Pi in some location, so that 

/\ = 01 = (23) 

Exploiting the above assumptions in Eq. (|22|) . we find that the Lagrangian £o(^) = K{t) of 
the ideal chain is given by: 

^oi^) = -r / ds{L — s) / du(p{t, s)(p{t,u) cos{(p{t, s) — (p{t,u)) (24) 
L Jo Jo 

What happens if, instead of polar coordinates, we choose cartesian coordinates in order 
to compute the continuous limit of the kinetic energy ([2])? With the help of the prescriptions 
given in Eqs. (fT2] - [T71) and related comments, it is easy to show that the Lagrangian of the 
ideal chain Co^disc = Kdisc becomes in cartesian coordinates: 

M rL 

Ut) = 7^ dsix\t,s) + y\t,s)) (25) 
ZL Jo 

Of course, the fields x(t, s) and y{t, s) are not independent, because they satisfy the rela- 
tionship: 

ix\t,s)f + {y'{t,s)f = l (26) 

where x' = || and y' = ||. Eq. ( |26l) is the continuous version of the constraints (jlj). As 
we see, the form of the kinetic energy is much simpler than that of its analogue in polar 
coordinates, but the price to be paid is the cumbersome presence of the constraints fl26l) . If 
li and (fii are constants, according to the assumption of Eq. (123!) . one should also add to 
Eqs. ( !25l) and ( l26l) the boundary conditions: 

x{t,0) = licosipi y(t,0) = lisimpi (27) 

In this way one end of the chain is fixed at the given point (x(t, 0), ?/(t, 0)) = 
(Zi cosy?!, Zi siny^i). It is possible to implement other boundary conditions. For instance, 
one could ask that the chain forms a closed loop: 

x{t,0)=x{t,L) y{t,0) = y{t,L) (28) 

To show that the Lagrangian in cartesian coordinates (l25l) and the Lagrangian in polar 
coordinates are equivalent, it is possible to perform in Eq. fl25|) the field transformations: 

x{t,s)= / du cos (p{t, u) + h cos ipi 
Jo 

y{t,s)= / du sin ip{t, u) + li sin (fi (29) 

JO 
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These are analogous to the discrete changes of variables of Eqs. ([6]) and ([7]) in the case in 
which one end of the chain is kept fixed. It is easy to check that, after the substitutions fl29|) 
in Eq. ( l25l) . one obtains exactly Eq. I^^ as desired. 

For future convenience we introduce also the vector notation: 

R{t,s) = {x{t,s),y{t,s)) (30) 

In this way we get for the functional Cq and the constraint (l26l) the more compact expres- 
sions: 

^oit) = ^l dsR\t,s) (31) 

and 

{R\t,s)f = l (32) 
Finally, Eqs. ( !27l) and ( !28l) become respectively: 

R(t, 0) = (/i cos ipi, h sin ipi) (33) 

R(t, 0) = R(t, L) (34) 

It is now easy to add the interactions. For example, let us suppose that the segments of 
the chain are immersed in an potential Vi(r) and that there are also internal interactions 
associated to a two-body potential V2(ri,r2). In this case, the Lagrangian Cq of Eq. (13T|) 
generalizes to: 

C = Co + Ci + C2 (35) 



where 



and 



£i = - dsVi{R{t, s)) (36) 
Jo 

C2 = - r ds^ r dS2Vint{R{t, Si), R(t, S2)) (37) 

Jo Jo 



III. DYNAMICS OF A CHAIN IMMERSED IN A THERMAL BATH 



In this Section it will be addressed the problem of describing the dynamics of a random 
chain subjected to thermodynamic fluctuations and immersed in an environment held at 
constant temperature T. 
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The strategy in order to treat this problem is to consider the discrete chain as a set of 
particles of mass m performing a random walk while subjected to the discrete constraints 
dl]), which we rewrite here for convenience as follows: 

|R„(t)-R„_i(t)|^ n = 2,...,iV (38) 

It is additionally required that at the instants t = ti,tf the n—th particle is located respec- 
tively at the initial point Ri^„ and at the final point R/,n for n = 1, . . . ,N. For sim plic ity, 



25| 



the interactions among the particles are switched off including hydrodynamic forces 

If one could ignore the constraints, the probability distribution \l/7v of the system of N 
particles would be: 

TV 

^iV = n Mtf - U, R/,n, Ri,n) (39) 

71=1 

where ipn is the probability distribution describing the free random walk of the n—th particle. 
As it is well known, ipn satisfies the Fokker-Planck-Smoluchowski equation 

= (40) 

D being the diffusion constant. Eq. (I40l) is completed by the boundary condition: 



^-71(0, R/,„, Ri,n) = 5(R/,n - Ri,n) (41) 

The solution of Eq. (140|) can be written up to an irrelevant normalization factor A in the 
form of a path integral: 

^n = Aj^^^^^^^^^ JR^it)e-k - (42) 

R-n(ti) = R-i,n 

Substituting Eq. fH21) in Eq. fl5^ . the probability distribution \I^jv becomes: 

= / ^ ^ dRiit) ■■■ I , , dRJt) exp I - y ?^dt\ (43) 

R-i(*i)=R-i,i R-jv(*i)=R-i,]v '"=1 ' 

Now we have to add to the above free random walks the constraints (|38l) . This will be 
done by inserting in the probability distribution of Eq. (H3l) a product of Dirac 5— functions 
which enforce exactly these constraints: 



R-i{ti)=R-i,i Rjv{t,)=R-i,jv 



R„(t) - R„_i(t)r 

(44) 
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In the above equation C denotes a normalization constant. Moreover, we have introduced 
the functional Ao^disc defined as follows: 

^^04^sc = E / (45) 

n=l "^^^ 

The distribution "^disc measures the probability that a chain starting from the initial con- 
figuration {Ri,i, . . . , Rj,Ar} of its segments ends up after a time tj — ti in the configuration 
{R/,1, . . . , R/^at}. We note that the diffusion constant D appearing in ( H5|) satisfies the 
relation: 

D = fxkT (46) 

where fi is the mobility of the particle, k is the Boltzmann constant and T is the temperature. 
This fact allows to rewrite the quantity Ao^isc in a form which reflects the analogy of the 
present problem with a quantum mechanical problem: 

1 A' 

Ik 



Ao,d.sc = ^ E I (47) 



n=l ' 

In the above equation we have put 

T = fim (48) 

The parameter r has the dimension of a time. Indeed, r is the relaxation time that charac- 
terizes the rate of the decay of the drift velocity of the particles composing the chain. The 
quantity Ao,disc looks now like a real action of a set of N quantum particles of mass m in 
complex time, with the constant factor 



K 



2kTT (49) 



replacing the Planck constant h. This is not a surprise, because the connections between 
quantum mechanics and Brownian motion are well known. Indeed, one may show that 
the uncertainties in the position and momentum of a Brownian particle are related to the 
constant 2mD = k 171 ]. 



At this point we are ready to take in the probability distribution (jH]) the continuous 
limit ( fT2l) . By introducing the rescaled variables s„ as we did in Section [Tll the probability 
distribution (144|) becomes: 

^^^sc = Cl[ ^^,.„,..^,„,c^R(M„)e-^4'^^-^^(*-)^^" 



n=2 



a 



13 



In the limit — > oo, a — > 0, Na = L we obtain from "^disc the probabihty distribution 
of the continuous chain: 

2 



where 



R(ti,s)=Ri(s) 



(9R(t, s) 



9s 



(51) 



This is the desired result. Formally, the normalization constant C may be written as a path 
integral over the initial and final configurations: 

2 



R(ti,s)=Ri(3) 



9R(t, s) 



ds 



1 (53) 



The model described by Eqs. ( 15TI - I52!) will be called here the generalized nonlinear sigma 
model due to its close resemblance to a two-dimensional nonlinear sigma model. The most 
striking difference is that the constraint = 1 of the nonlinear sigma model has been 
replaced by the condition fl5^ . which contains the derivatives of the bond vectors R and it 
is thus nonholonomic. 

Before concluding this Section, we would like to complete the derivation of the probability 
distribution of Eqs. fl5T] - [52l) by discussing the continuous limit of the relaxation time r. This 
parameter has been defined in Eq. ( I48l) as the product of the mobility fi with the mass m 
of the beads. When the distance a between the beads goes to zero, m goes to zero as well 
according to Eq. (IT^ . On the other side, with decreasing values of a, two contiguous beads 
will become closer and closer until they will eventually merge one into another. To avoid 
this unphysical situation, one should add to the continuous limit f|T2|) the requirement that 
the dimensions of the beads vanish together with a. Supposing for instance that the beads 
are circles of radius p, for dimensional reasons one is lead to put p = ca, where c is a 
dimensionless proportionality factor. As it is intuitive, when the size of a bead decreases, 
its mobility p increases. The increasing of p compensates the vanishing of m, so that the 
product r = pm remains finite. This fact can be verified rigorously in three dimensions 
using the well known Stokes formula of the mobility. 
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IV. THE CLASSICAL SOLUTIONS OF THE GENERALIZED NONLINEAR 

SIGMA MODEL 

To get rid of the delta function appearing in the formulation of the generalized sigma 
model of Eqs. ( l5Ti) and ( 152|) . it will be convenient to introduce a Lagrange multiplier A = 
A(t, s) and to use the Fourier representation of the Dirac delta function. Moreover, we 
add a coupling of the bond vectors R(t, s) with an external source J{t, s). In this way the 
generating functional "^[J] of the correlation functions of the bond vectors may be written 
in the form: 



^[J] = J VRVXexp I- J^^ dtj^ ds 



^ R2 + a(R"-l) + -J-R 



2Lk ^ ^ K 



(54) 



Let us note that, for the sake of generality, in Eq. ( l54l) no boundary conditions for the 
relevant fields have been specified. As it stands, Eq. f l5^ could be applied both to open or 
closed chains. Moreover, in the case of open chains both possibilities of free or fixed end 
points are allowed. It turns out that the degree of complexity of the computation of \Ef[J] 
strongly depends on the choice of boundary conditions. 

Let's now derive the solutions Rc; and Ac/ of the classical equations of motion associated 
with the generating functional '^[J] of Eq. ( l54l) : 

- J (55) 

^^ + A^ - (56) 

ds ds ds"^ 

1 (57) 



dR ' 



ds 

It is easy to see that, due to the constraint ( 1571) . the current J must depend only on the 
variable t, i. e. J{t,s) = J{t). Always for the same reason, it is possible to check that 
Eq. ( l56l) is inconsistent unless A^ = const 7^ or A^ = 0. Depending on the fact that A^ is 
zero or not, one finds that the relevant solutions of Eqs. fl55] - [57j) may be divided into two 
groups, which we call here solutions of type A and of type B. The solutions of type A are 
characterized by the condition = 0. They are displayed in Table [H Besides the classical 
solutions of type A, there are also the solutions of type B listed in Table [Tll 

We would like to stress the fact that, if there are no external currents, both type A 
and B of classical solutions are static, i. e. they do not depend on time apart from the 
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Type A solutions 



Ad = (58) 
Rd(t,s) = Ro,A + Vt + Ri,A(t) + R2,A(s) (59) 



where Rq a and V are constant vectors, 



Ki^Ait) = r dt'G{t,t')3{t') (60) 



and 

rs 



R2a(s)= / du{cos ip{u), sin ip{u)) (61) 
Jo 

Here (p{u) is an arbitrary function of u, while in Eq. (|60p G{t, t') denotes the Green function 
which solves the differential equation: 

,62, 

Type A solutions admit closed chain configurations. In that case, the functions 99(5) must satisfy 
the additional periodicity condition: 

ip{s + L)= 99(5) (63) 



TABLE I: Solutions of type A of the classical equations of motion (j55H57p . 

rigid translations of the whole chain with constant velocity V. This result is confirmed if 
one studies the equations of motion corresponding to the Lagrangian 0241) obtained using 
polar coordinates. The only allowed classical solutions for the field (p{t, s) are in fact time 
independent. 

V. COMPUTATION OF THE GENERATING FUNCTIONAL ^[J] IN THE 
SEMI-CLASSICAL APPROXIMATION 

The exact computation of '^[J] is a formidable problem despite the simplicity of the action 
of the generalized nonlinear sigma model. One of the main difficulty is the presence of the 
nonholonomic constraint |R'P = 1 in the path integral (ISTll . In principle, this cumbersome 
condition may be easily eliminated by introducing a scalar field (p{t, s) and performing the 
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Type B solutions 



Xci = const / (64) 
nd{t,s) = Ro,B + Vt + Ri,B(t) + R2,ij(s) (65) 



where Rq b and V are constants vectors, 



lii^Bit) = r dt'G{t,t')J{t') (66) 



and 

R2,ij(s) = (0,s) (67) 

In the absence of the external source J(t), this kind of solutions corresponds to a configuration in 
which the chain is stretched along the y— axis with one end in the point Ro,_b and the other end 
in the point Ro,b + (0, L). No closed loop configuration is allowed. Of course, the stretched chain 
may be oriented in a different way by means of a rotation. 



TABLE II: Solutions of type B of the classical equations of motion (|55H57j) . 
formal substitutions of Eq. (12^ 

R(t, s) = / du{cos (p{t,u), sin (p{t,u)) (68) 
Jo 

Here we have assumed for simplicity that the chain has one fixed end in the origin of the 
coordinates, so that R(t, 0) = (0,0) in agreement with Eq. ( l68l) . To show that after the 
field transformation fl68|) the constraint disappears from the path integral flSTl) . we use the 
following relation which will be proved in Appendix A for a generic functional /(R(t, s)): 

J PR(t,s)/(R(t,s))5(|Rf - 1) J Vip{t,s)f{£ du{cosip{t,u),smip{t,u))) (69) 



where is an irrelevant constant. In our particular case, in which: 



/(R(t,s)) = exp 
one obtains from Eq. (!69l) : 



dt I ds^R'^ 
Jo 2nL 



(70) 



^=U Jv^{t,s)e'^^^''''''' (71) 
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with £o being the Lagrangian defined in Eq. (12^ . As we see, the Dirac delta function with 
the constraint is no longer present in the path integral (171]) . but the Lagrangian of Eq. flM|) 
is both nonlocal and nonlinear. 

In the following, we will stick to cartesian coordinates limiting ourselves to study small 
gaussian fiuctuations of the field R(t, s) around the classical solutions derived in the previous 
Section. To this purpose, in Eq. fl54|) we split both fields R(t, s) and X(t,s) into classical 
contributions HchKi and statistical corrections 5R, 5A: 

K{t,s) = Rc/(t,s) + K^5R(t,s) (72) 
X{t,s) = \ci + KhX{t,s) (73) 

Moreover, it will also be convenient to split the external source J in an analogous way: 

Jit,s) = Jci{t)+Khjit,s) (74) 

where 3d(t) denotes the current depending only on the time t appearing in Eqs. (!60l) and 
fl66|) . Due to the fact that the Lagrange multiplier A is just an auxiliary field, it is possible 
to choose for its variation 6X trivial boundary conditions at the initial and final instants: 

6X{ti,s) =6X{tf,s) = (75) 

The boundary conditions of 5R(t, s) will be fixed later. 
At this point, we expand the action 

M /„,2 .\ 1 



l-tf rL 

A 



appearing in the path integral ( ]54l) with respect to the quantities 5R, 6X and 6J. Since the 
latter are supposed to be small corrections of the dominating classical solutions, we stop the 
expansion at the second order: 

A = A,i + 6A^^^ + (77) 

At the zeroth order we have: 

T rtf rL r M . ^ 1 

(78) 



1 /■*/ 

Ad = — I dt j ds 



2L 



'^^d + JcZ ■ Rd 



K Jti JO 

This is just the action A in which the fields have been replaced by their classical configu- 
rations, which may be either of type A or of type B. In both cases, the term Aci(R'^ — 1), 
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which in principle should be present in Eq. ( !78|) . has been omitted because Eq. (l57|l forces 
it to vanish identically. Let's now compute the first order contribution SA^^''. Usually, first 
order contributions vanish after exploiting the classical equations of motion. In our case 
this is in general not true. The reason is that, due to the nontrivial boundary conditions 
satisfied by R(t, s), nonzero boundary terms may still appear in SA^^\ Despite this fact, it 
is possible to show that 5A^^^ vanishes at least in the following two situations: 

1. Closed chains satisfying the boundary conditions 

2. Open chains in which both ends are fixed, so that besides the condition fl33|l also the 
following one is valid: R(t, L) = (xq,/,, ?/o,l); where Xo,l and |/o,l are constants. 

Assuming that one of the above two conditions is verified, it is possible to put: 

SA^^^ = (79) 

Thus, we are left only with the computation of the second order corrections SA^'^\ After 
simple calculations, one finds: 

5^(2) = f'^ dt ds [^5R ■ 5R + inX^SR'Y + 2m(R'^^ ■ SR')5X + 53 ■ 5r] (80) 
Jti Jo i2L J 

Putting together all the above results, it is possible to conclude that, within the present 

gaussian approximation, the expression of "^[J] reduces to: 

^[J] = e-^^'ZiJ] (81) 

where 

Z[J] = J V6RV6Xe-^^^^'^ (82) 

and Ad and 5A^'^^ are respectively given in Eqs. (ITS!) and (IHOj) . 

From this point on we will consider only chain configurations which are closed, so that 
both classical configurations Rd and their statistical corrections 6R must satisfy the bound- 
ary conditions in s of Eq. ( IMl) . In the Lagrange multiplier sector, besides the trivial bound- 
ary conditions in time of Eq. (1751) . we require also the following ones with respect to the 
variable s: 

6X{t,0) = 6X{t,L) (83) 

The case of closed chains is particularly interesting because under this assumption the path 
integral appearing in the right hand side of Eq. (I82l) may be rewritten in such a way that it 
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closely resembles the generating functional of a field theory of a one- dimensional system at 
finite temperature. In this field theory the coordinate s plays the role of time while the real 
time t becomes the spatial coordinate of the one-dimensional space. We still need to specify 
the boundary conditions with respect to the time for the fields SH. We fix them in such 
a way that the gaussian path integration over these fields in the generating functional of 
Eq. (182|) will be as simple as possible. To this purpose, reasonable choices are the following: 

Dirichlet-Dirichlet boundary conditions 



SR{ti, s) = 



5R{tf,s) = 



(84) 



Dirichlet-Neumann boundary conditions 

SK{ti, s) = 

Neumann-Dirichlet boundary conditions 

d6R{t, s) 



dt 



t=ti 



Neumann-Neumann boundary conditions 

d6R{t, s) 



dt 



t=ti 



(9R(t, s) 



dt 



t=tf 



6R{tf, s) =0 



dSR{t, s) 



dt 



t=tf 



(85) 



(86) 



(87) 



At this point we are ready to perform the Gaussian integrations over the fields 6R in the 
generating functional Z[J] given in Eq. fl52]) . After some integrations by parts, which do 
not produce boundary terms thanks to the boundary conditions and one finds: 



Z[J] =CiJ I?5Ae-^(^^) 



(88) 



with 



S{5\) = I ds dtdt'Git, t' 
2 Jo Jti 



— 2iK6X{t, s] 



-4K^5A(t,s 
dR,i{s) d63{t',s 



6X{t', s 



ORciisY 



ds 



ds 



- 2m5\{t\ s) 



ds ds"^ \ ds 
dR,iis) dSJit,s) , 



ds 



ds 



+ 5J{t,s) ■ 63{t',s) 



(89) 
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and 



Ci = J P5Re'^Jo (90) 
The symbol Git, t') denotes the propagator ([52]) computed taking into account one of the 



boundary conditions defined in Eqs. flMH Hl 
Dirichlet-Dirichlet boundary conditions 



18|. We have thus four possibihties: 



G(t,t') = -^e{t' -t){t~ u)^^^—^ - - t'){t' - u)'^-^ (91) 



Dirichlet-Neumann boundary conditions 



G{t,t') = -^{u - mt - 1') - - mt' - 1) (92) 



Neumann-Dirichlet boundary conditions 



G{t, t') = ^{tf- mt' -t) + ^{tf- t)e{t - t') (93) 
Neumann-Neumann boundary conditions 

In Eqs. (!9T] - [9^ the function 9{t) is the usual ^—function of Heaviside. Let us note that the 
function G{t, t') of Eq. flMl) is actually a pseudo Green function, which satisfies the equation: 

Md^Git^f) „ 1 

instead of Eq. fl62l) . This is due to the fact that, in the case of Neumann-Neumann boundary 
conditions, one should project out the constant solution of the eigenvalue equation associated 
to the operator 

We remark that in the action (l89l) the classical conformations appear only in the deriva- 
tives which coincide with the derivatives of the vectors R2,yi("S) defined in Eq. (I6T!) . In 
components: 

= (cosv2(s),sinv9(s)) (96) 

From Eq. 0961) it is clear that the vector has the meaning of the unitary vector which 
is tangent to the classical trajectories. It is therefore convenient to decompose all vectors 
appearing in the action S{5\) in components which are normal or tangent to Rc/(s). After 
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some algebra, one obtains in this way an expression of 'Z\J\ in which now S[bX) takes the 
simphfied form: 



1 



S{b\) = I ds dtdt'G{t, t') 
2 JO Ju 



(97) 



+ 2iK6X{t, s)6JT{t', s) + 2iK6X{t', s)6JT{t, s) + 63{t, s) ■ 63{t', s) 



Here we have introduced the tangential component SJxit, s) of the current 63{t, s): 



SJrit, s) = dsRciis) ■ ds6J{t, s) 



At this point we are ready to ehminate the auxihary field 6X from the functional of Eq. (1881) . 
where now the action S{6\) is defined in Eq. ( 1971) . To this purpose, one needs to perform a 
gaussian integration, which produces the result: 



(98) 



(99) 



In the above equation K{s, s') denotes the Green function which satisfies the relation: 



K{s, s') = —S{s, s') 



while 



ds f^f dtdt'G{t,t') 2K^5\{t,s)(^j^-{^'{s)f'^&\{t' ,s) 



(100) 



(101) 



Putting all together we obtain the expression of the generating functional '^[J] in its final 
form: 



$[J] = e-^^'CiCaexp 

1 rL 
X exp -- / 

2 Jo 



\ ds dtdt'Git, t')S3{t, s) ■ SJ{t', s) 
L Jo Jti 

dsds' dtdt'G{t, t')K{s, s')6JT{t', s'))6JT{t, s) 
Ju 



(102) 



The right hand side of Eq. (11021) displays the asymmetry in the propagation of transverse 
and longitudinal modes. 

The differential equation satisfied by the Green function K{s, s') of Eq. fllOOl) cannot be 
solved analytically for any given function ip{s). Here we discuss just the case in which the 
background classical solution corresponds to a chain with the configuration of a circle, i. e.: 
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The radius of the circle is so that the total length of the chain is L as desired. Comparing 
Eq. (11031) with Eq. (!6T|) . it is clear that for this conformation ip{s) = ^s. Substituting this 
expression of ip{s) in Eq. (llOOp . it turns out that the Green function K{s,s') satisfies the 
relation: 



An' 



K{s, s') = -6{s-s') 



The solution of the above equation corresponding to the boundary conditions (J83ll is: 



Kis,s') 



sinh 



^(L-s')l sinh^s 



^ sinh 27r 



e{s' -s) + 



sinh 



f(L 



^sinh27r 



(104) 



sinh ^s' 

■^e{s - s') (105) 



In the limit J(t, s) = 0, we obtain from the generating functional of Eq. (I102p the expres- 
sion of the probability distribution \1/ of Eqs. fl5T| - [52l) in the semiclassical approximation: 

(106) 



Remembering the respective definitions of the constants Ci and C2 of Eqs. (1901) and (IIOII) . 
together with the form flTHj) of the classical action, may be explicitly written as follows: 



\1/ = exp 
X / D^Aexp 



— / dt I ds 

K Jti Jo 



M. 



R 



V6H exp 



2L Jti Jo 



ds (5R)' 



*^ dtdt' r dsG{t, t') \2ti^5\{t, s) (£j - {v'{s)f) 5X{t', 

ti JO 



(107) 



VI. PHYSICAL INTERPRETATION OF THE OBTAINED RESULTS 



The classical equations of motion (!55] - [57I) admit only static solutions in which the con- 
formation of the chain is fixed. Only rigid translations of the chain as a whole with constant 
velocity V are allowed. Apart from these rigid translations, the time dependence of the 
classical solutions Hci(t,s) in Eqs. fISUl) and flB5]) is just an artifact of the presence of the 
classical current J{t) = Jci{t) which is fictitious and may thus be set to zero without any 
loss of generality. In the following discussion, it will be assumed that both 3d{t) and V are 
zero. The absence of any relevant dynamics in the classical solutions is somewhat surprising. 
We can only suggest that this absence could be related to the fact that, after performing 
the continuous limit, the fourth term appearing in the left hand side of Eq. flTT]) vanishes 
identically. This term is important for the chain dynamics since it contains second time 
derivatives of the angles ip{t,s). In some sense, the continuous chain is simpler that its 
discrete counterpart, in which this term is present. 
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Less trivial is the treatment of the fluctuations of the chain at constant temperature T. 
The nonhnear sigma model given in Eq. fl5^ or, alternatively, its formulation without the 
Lagrange multiplier of Eqs. ( I5TI - I521) . describe the fluctuations of a chain of segments of 
constant length a in the continuous limit ( |T2|) . During its motion in the time interval 
the chain spans a two-dimensional surface which, in the case of a closed conformation, has 
the topology of a cylinder in the x, t space, see Fig. [2l If the chain is open, instead, the 




FIG. 2: During its motion, a closed chain spans in the x,y,t a surface which has the topology of 
a cylinder. 

topology of the cylinder should be replaced with that of a strip. 

We recall that in deriving the GNLSM of Eqs. ([5IHS2D the contribution of the hydrody- 
namic interactions has been neglected. This limits the validity of this model to the following 
cases: 

a) The viscosity of the fluid is large, so that the motion of the particles composing the chain 

is slow. Thus the velocity fleld created by each particle is too weak to influence the 
motion of other particles. 

b) The temperature is low, so that once again the motion of the chain is slow. 

c) The conformation of the chain is relatively straight because there is some energy cost 

when the chain is being bent. To this purpose, however, one should introduce the 
stiffness at the joints between the segments. 

On the other side, the semiclassical approximation used in order to derive the generating 
functional "^[J] of Eq. fll07p is valid in the case in which the parameter k, deflned in Eq. fl49p is 
small. This parameter depends essentially on the temperature T and on the relaxation time 
r. Since r is inversely proportional to the viscosity in the limit of low Reynolds number, it is 
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reasonable to assume that the semiclassical approach can be apphed to a cold isolated chain 
or to an isolated chain in a very viscous solution. These situations correspond respectively 
to the points b) and a) mentioned above. 

Both the generating functional "^[J] of Eq. f ll02p and the probability distribution of 
Eq. fll07p have been computed in the case of closed chains, whose conformations are subjected 
to the boundary conditions (jHHSTj). The physical meaning of these boundary conditions 
may be summarized as follows. 

Dirichlet-Dirichlet boundary conditions: In this case Ri(s) = R/(s) = Rci(s), where Rc/(s) 
is a given static solution of the classical equations of motion. The probability that, 
starting from the conformation Ri(s), the fluctuating chain ends up at the time t/ in 
the same conformation, is proportional up to a normalization constant to the prob- 
ability distribution \l/ of Eq. fllOTp . in which the Green function G{t,t') is that of 

Eq. dnn). 

Dirichlet-Neumann boundary conditions: In this case the probability distribution \1/ of 
Eq. fll07p is proportional to the probability that a closed chain starting from a classi- 
cal static conformation Rj(s) = Rc/(s) at the time ti ends up at the instant tj in an 
arbitrary conformation characterized by the fact that the velocities of each segment 
composing the chain is zero. \1/ must be computed choosing the Green function G{t, t') 
defined in Eq. ((921). 

Neumann-Dirichlet: Here the segments of the chain have zero velocity at the beginning, but 
the conformation of the chain is otherwise arbitrary. At the time tf the chain is found 
in a given static classical conformation, i. e. R/(s) = Rci(s). The probability for this 
to happen is obtained after substituting in Eq. (I107P the Green function G{t, t') of 
Eq. (I93D. 

Neumann-Neumann: This is the situation in which the conformation of the chain at the 
initial and final times ti and tf are not specified, but the velocities of all the segments 
composing the chains must be zero. The relevant Green function t') to be inserted 
in the probability distribution \1/ is that of Eq. (IMjl . 

By choosing Neumann-Dirichlet boundary conditions one may check for instance if, starting 
from any static conformation, there is a particular conformation in which it is very likely 
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that the chain will be found after a certain time tf — ti. The stability of a given conformation 
with respect to the thermodynamic fluctuations which attempt to reshape the chain can be 
tested by choosing Dirichlet-Dirichlet boundary conditions. In principle, it is also possible 
to study other types of boundary conditions than those considered here, provided they do 
not give rise to unwanted boundary terms in the action of the generating functional '^[J]- 

To conclude this Section, it is interesting to see how the original constraint (l32l) is realized 
in the semiclassical approximation. In Eq. flHTl) . which gives the second order corrections to 
the classical action Acu there are two Lagrange multipliers, 5\ and Ac;. The most relevant 
condition is that imposed by 5X: 

■ (5R' = (108) 

Let us note that the above relation is at the origin of the asymmetry in the generating 
functional (11021) between modes which are tangent or normal to the classical background 
conformation R^;. Eq. (11081) is just the approximated version of the full constraint fl32l) . As 
a matter of fact, remembering the splitting into classical solutions and statistical corrections 
of Eqs. ( 1721) and ( 1731) . we may rewrite Eq. ( l32l) as follows: 

(R')^, + 2R;, ■ 5R' + {5^f -1 = (109) 

Due to the fact that (R'c;)^ = 1 and neglecting the second order term ((5R')^, we obtain from 
Eq. (UnHl): 

2R;,,-5R' = (110) 

which coincides exactly with Eq. (11081) . In the case of solutions of type B there is an 
additional constraint, which is associated to the nonzero constant Lagrange multiplier Xci. 
This constraint requires that the average over the time t and over the chain length s of the 
quadratic term ((5R')^ is zero: 

/ T rdt -(5R')' = (111) 

Jti tf — ti Jo L 

In the solutions of type A the above condition is not present because in that case Ac/ = 0. 
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VII. THE EQUILIBRIUM LIMIT OF THE GNLSM AND ITS CONNECTION 

WITH THE ROUSE MODEL 

First of all, we study the equilibrium limit of the GNLSM. We use the formulation of the 
model given in Eqs. (15T1452I) . For simplicity, we set ti = 0. Thus 

^ = J I?R(t,s)e-^i'o''^*/o^'^^^'(*'^)5(|R'(t,s)|' - l) (112) 

At this point, we rescale the time variable by putting a = ^, so that the above equation 
becomes: 

^ = y PR(tja,s)e'^^« '^"^o "''^^^ 5 (|R'(t/(T,s)|' - l) (113) 
In the equilibrium limit ^ oo we obtain: 

V^e, = / VR,,is)5 (|R'e,(s)f - l) (114) 

where Reg(s) = R(cx3,s). Eq. flll4p is exactly what one should expect in the case of the 
statistical mechanics of a discrete chain subjected to the constraints: 

|R„-R„_i|' = (115) 
The discrete probability function of the conformation of such a chain is given by: 

N N /IP} _ pj |2 \ 

^e,4^sc = n / n M ' f"' " 1 (HG) 

n=l 71=2 \ 0, J 

In the continuous limit this becomes exactly the distribution of Eq. 01141) . This result is in 
agreement with the analogous probability function given in [2]. 

While the purpose of this work is to provide a path integral formulation of the dynamics 
of a freely jointed chain without having in mind concrete applications to polymer physics, it 
is interesting to explore possible connections between the GNLSM of Eqs. ( 15T] - [52|) and the 
Rouse model. A direct attempt to put the Rouse model in the path integral form via the 
Martin-Siggia-Rose formalism leads to a probability distribution for the Rouse chain which 
differs profoundly from the GNLSM obtained in this work. Indeed, let us start from the 
Langevin equation: 

,9R (92R ^ , ^ 

Here we have used instead of the arc-length s the dimensionless variable n defined as follows: 
Son = s. Moreover, f = f (t, n) is a stochastic force with a Gaussian distribution of width a 
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given by: e Jo Jo "2". C cind ^ are constant parameters which will be specified later. 
After the application of the Martin-Siggia-Rose method, one finds the Rouse probability 
distribution: 



Rouse 



VRe 



2/' aRN^^ 

dt ) 



2 ( a^R- A ' 



(118) 



In principle, in the exponent of the above equation there should be the additional term 



a Jo 



(119) 



However, due to the fact that # ■ f# = # f ^ ■ - ^ ■ and remembering the 



dn 



dn \ dn dt J dn dtdn 



identity ^ 



dn dtdn 



2 dt 



(my 



\ dn ) 



it is easy to realize that I amounts to total derivative 
terms, which can be neglected. 

Coming back to Eq. (11181) . we see that, while the GNLSM is nonlinear and contains just 
second derivatives of the bond vector R, the Rouse model is linear and contains derivatives of 
R up to the fourth order. As an upshot, while it is possible to investigate the Rouse model 
by decomposing R(t, s) into normal coordinates as explained in [l], that kind of Fourier 
analysis cannot be easily applied to the nonlinear GNLSM. We show at this point that, 
indeed, the two models are quite different and that it is not possible starting from one of 
them to recover the probability function of the other and viceversa, because they correspond 
to regimes which do not overlap. To this purpose, instead of the constraint R'^ = 1 of the 



GNLSM, we introduce the more general condition 



7^ 



So "'-so 

where Sq is a new length scale such that 



1 /"^o 

- / ds'A{s') (R(t, s + s') - R(t, s)) 

5n J —sn 



— a 



a <^ Sq <^ L 



(120) 



(121) 



;i22) 



and A{s') is a function of s' normalized in a such a way that: 

- / A{s')ds' = 1 

So J-so 

Let us note that a is the smallest length at our disposal: Any segment of the chain of contour 
length shorter than a may be regarded as rigid. At this point, following the same strategy 
used in Section llllt we build the new distribution function 



(123) 
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with c = 4^2VL' index int means that the distribution probabihty \l/j„f describes a model 
which, as we will see, interpolates between the GNLSM and the Rouse model. 

We remark that the insertion of the 5-function S{TZ) in the path integral fll23p has a 
double meaning. On one side, it may be seen as a condition on the length of the averaged 
vector: 

S = — / ds'A{s'){R{t,s + s')~R{t,s)) (124) 

So J-so 

In the above equation the distance between points of the chain has been averaged over arc- 
segments of length 2so. On the other side, the introduction of the 5-function S{TZ) may also 
be related to the internal forces among the beads composing the chain, which appear due to 
the presence of constraints. For example, in the case of the GNLSM the presence of these 
forces is evident in the formulation of Eq. ([SID, which the free action is corrected by the 
addition of the interacting term A(R'^ — 1). 

It is easy to realize that the GNLSM is a special case of the model described by Eq. (11231) . 
We have just to remember that in the GNLSM the motion of the chain is observed at the 
smallest available scale of distances, i. e. the segment length a. Thus, we choose the form 
of the function A{s') as follows: 

A{s') = so6{s' - a) (125) 

As a consequence, the constraint (I120p becomes |R(t, s + a) — R(t, s)f — = 0. Dividing 
both members of the above equation by and supposing that a is very small, we get up to 
higher order terms in a the relation: 

|R'|^-1 = (126) 

In the limit a = 0, this is exactly the condition which has been imposed to the chain in the 
GNLSM, see Eq. ([321). Using the property of the 5-function ^(a^dR'p - 1)) = ^5(|R'|2-1), 
it is possible to check that also the probability distribution \E'j„t becomes that of the GNLSM. 

To obtain the Rouse model from the interpolating probability distribution "^int, we need 
first of all to decrease the resolution with which the segments of the chain are observed. 
Accordingly, we require that the function A{s') appearing in the constraint (11201) is constant 
over the whole interval [—Sq, Sq]: 

A{s') = 1 (127) 

In this way, the finest details of the chain are not taken into account, because the chain 
conformations are averaged over the scale of distance 2so, which is by hypothesis much 
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larger than the smallest available scale a. To pass to the Rouse model, we have also to 
restrict ourselves to the long time-scale behavior of the chain. This is achieved by assuming 
following Section 4.1, p. 93, that, for long-time scales, R(t, s) varies slowly with s. This 
hypothesis allows to stop the expansion of Hit, s + s') with respect to s' at the first few 
orders: 

R(t, s + s')= R(t, s) + R'(t, s)s' + 5l^|lfls'2 ^ _ _ _ (^28) 

Substituting the above truncated expansion in Eq. (11201) and performing the trivial integra- 
tions over s' we obtain the condition: 



= (129) 



62 

Plugging in the above constraint in Eq. fll23p , we get the following approximated expression 



of ^ 



int 



^int ~ / VRe-'Io' '''Io'''^'6{\K"\hl - «') (130) 

where the factor 6^ has been absorbed by a rescaling of the length sq- At this point we use 
the fact that, apart from an irrelevant infinite constant, the functional 5-function present in 
Eq. (11301) may be simplified as follows: 

5(|R"|'s^ - a') = 5(|R'>o - «) (131) 

A proof of this identity, which is valid up to an irrelevant constant, will be given in Ap- 
pendix B. Exploiting Eq. (11311) . the probability distribution "^int of Eq. (I130p becomes: 

~ / PRe-^/o''^*/o"'^^i^^5(|R>2 _ ^) (132) 

We may still simplify the above equation by applying the following slightly modified version 
of the gaussian approximation of the (5-function: 



6{\R \st,-a) J VXe -"^ -"^ V' '° J e ' ^ ) (133) 

where we have supposed that the parameter v is very large while the ratio (3/v is very small. 
Clearly, the usual Fourier representation of the functional Dirac ^-function is recovered in 
the limit v 00 and (3/v — > 0. Up to now (3 is an arbitrary parameter. We choose it in a 
such a way that 

^ = a (134) 
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This choice is compatible with our requirement for /3, since a is the smallest scale of lengths 
at our disposal, so that P/u is a very small quantity. Using Eq. (11341) in order to eliminate 



u from Eq. (I133p . we obtain the relation: 

5(|R'>^ _ a) ~ I VXe-'Io'^'Io"^^'\^"\4e-i:/'^^I>^^^' (135) 

After performing the gaussian integration over A in Eq. ( 11350 . we obtain 

d)(|R |sq — a)~e Jo Jo 2ai i o (136) 

We may now plug in the above expression of the 5-function in the distribution probability 
of Eq. ffT^ . The result is: 

^,„t ~ / VRe ^0 Jo [ 2.1 I oj ^^3-^^) 

This approximated probability distribution has the same structure of that coming from the 
Rouse model given in Eq. (IllSp . 

To make the comparison with the Rouse model more explicit, we perform in Eq. ( 1137P 
the substitution nsQ = s: 



T / ^T-> Jo Jo o 1 2a a„2 r 

"^int ~ J VRe I J 



(138) 



Let's now identify the coefficients appearing in Eq. (I138P with those of Eq. (11180 . We recall 
the fact that in the case of the Rouse model: 

C^i a^i^ (139) 

On the other side, the parameter c is in the exponent of Eq. (I138P may be written as follows: 
c = It is now easy to verify that the probability function Rouse of Eq. (I118p and that 
of Eq. (I138P coincide if we make the following choice for /3: /3 = |^/cT. 



VIII. CONCLUSIONS 



This work may be considered as the ideal continuation of the seminal paper of Edwards 

n 

and Goodyear of Ref. [2], in which the problem of a chain subjected to the constraints ( 138|1 
has been investigated using an approach based on the Langevin equation. With respect 
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to Ref. [2^, our approach based on the Fokker-Planck-Smoluchowski equation provides a 
path integral and field theoretical formulation of the dynamics of a freely jointed chain in 
the continuous limit. The GNLSM obtained here makes possible the application of field 
theoretical techniques to the study of the fiuctuations of a freely jointed chain. As an 
example, we have derived in the semiclassical approximation the probability function of the 
chain and the associated generating functional. The approximation used in the computation 
of the generating functional is valid for instance in the cases of a cold isolated chain or of a 
chain fiuctuating in a very viscous medium. 

Most of the results obtained in this paper have been discussed in the previous two Sections, 
so that we provide only a short summary of them: 

1. Derivation of the GNLSM, which provides a path integral formalism to the freely 
jointed chain. 

2. Computation of the partition function and of the generating functional of the GNLSM 
in the semiclassical approximation. 

3. The behavior of a chain at scales of length and time which are very long has been 
compared with the behavior at short scales of length and time in Section VII. It is 
shown in this way that it is not possible to compare directly the Rouse model and 
models describing a freely jointed chain like the GNLSM, because the regimes and 
the assumptions of these two models do not overlap. In particular, the Rouse model 
considers only the long-time behavior of the chain and long scales of distances, while 
in the case of the freely jointed chain the short-time behavior is taken into account 
and the chain is observed at a short scale of distance. 

4. A chain model which encompasses both the regimes of the Rouse model and of the 
GNLSM has been proposed in Eq. ffT23|l . 

5. The equilibrium limit of the GNLSM has been recovered. It gives the expected result 
in agreement with Ref. Q]. 

6. Last but not last, the dynamics of a random chain has been investigated also from the 
classical point of view. The equivalence of the expressions of the classical Lagrangian 
of the chain computed starting from cartesian and polar coordinates has been verified. 
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To conclude, we would like to mention the problems which are still open and possible 
further developments of our work. For simplicity, we have restricted our analysis to a two- 
dimensional chain. However, its extension to any dimensions is not difficult. Preliminary 
work in three dimensions can be found in Ref. lol. It turns out that more dimensions 
allow the possibility of performing the continuous limit in different ways, so that one could 
end up with a flexible chain or with a rigid chain which privileges only certain directions 
along a fixed axis. There should be also no problem in switching on the interactions among 
the beads composing the chain. To this purpose, one may use the path integral methods 
applied to stochastic differential equations explained in Ref. js]. Only the inclusion of the 
hydrodynamic interactions requires still some work. It is not simple to provide for these 
interactions a lagrangian based formulation as that developed here for a continuous chain. 
However, hydrodynamic interactions have been already implemented in the path integral 
formalism in Ref. 20] u sing the Martin-Siggia-Rose formalism. Work is in progress in order 



to extend the results of 201] to the freely jointed chain discussed in this work. Another open 
question is how the chain behaves in the short time regime when it is stretching under the 
action of a force. This could be interesting in the biophysics of DNA 2l|]. Finally, work is 
in progress in order to linearize the GNLSM applying the same approach used in the case of 
the standard nonlinear sigma model. In this way it would be possible to study the GNLSM 



as the strong coupling limit of its linear version with the help of the techniques of Refs. 



22]. 
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APPENDIX A: PROOF OF EQ. (f69]> 



In this Appendix we wish to prove Eq. (169!) . To this purpose, we start from the path 
integrah 

I = J Vx{t, s)Vy{t, s)f{x{t, s), y{t, s))6 [{d^xf + {dsvf - l) (Al) 
Upon the transformation: 



Xs{t,s) = dsx{t,s) 



(A2) 



we obtain: 



/ = J Vxsit,s)T>y,{t,s) (det ^ds^ f (^J^ duxuit,u), duyuit,u)^ S (^x^ + - 1^ 
where we have made use of the fact that: 

Now it is possible to ehminate the variable ys by performing in / the substitution: 



(A3) 



det ( 1 det (-^ 
\dxsl \dys 



(A4) 



X = xl + yl-l 



(A5) 



As a consequence, the path integration over ys appearing in ( ]A3I) may be replaced by a path 
integration over the new variable x'- 



J Vys5 (xl + yl-l^ f l^j^ Xudu, yudu 



Vx5{x) det 



f (^J^ Xudu, yuduj 



(A6) 



In the above equation the determinant det ^ is the functional determinant giving the 
Jacobian of the transformation ( 1A5I) . so that: 



det 



det 



2a/1 + X-x2 



(A7) 



Applying Eqs. (Ell) and flMl) in flXTj) we get: 



Vxs det \d/ det 27l - / / u)du, ± Jl- xl{t, u) (A8) 
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Finally, we perform in Eq. (lASP the substitution Xg = cos (f. In this way I may be rewritten 
as a path integral over ip: 

5x. 



Dip det 



Noting that det 



Sip 

SXs 



Sip 



det ^ 9g det ^ {2sm(p) f (^J cos (p{t,u)du, J sm(p(t,u)diLj (A9) 
= — sin (f we obtain, apart from an irrelevant constant A/" = 



det ^ |9g I det ^2, the final result: 

I = J\f J Vipf (^J cos (p{t,u)du, J sin Lp(t,u)du 



(AlO) 



APPENDIX B: A FUNCTIONAL IDENTITY 



To prove the identity fll3ip . we write the functional 5— function in the right hand side of 
Eq. (11311) with the help of its Fourier representation: 



5 (iR'fst -«')=/ V~Xe-'Io'^^I>^~^(\^"\'<-') 



(Bl) 



Since \R."\h^ - = (|R"|s^ - a){\R"\sl + a), Eq. 1^ becomes: 

6 (|R'f ^^e-*/o' <it/>.A(lR"|.g-a)(|R"lsg+a) (32) 

Due to the fact that |R"|so + a is always different from zero, one may perform the change 
of variables: 

A=~X{\R"\sl + a) (B3) 



Applying the substitution (IB3p in (IB2p we obtain: 

S (|R'f - a') = det-^(2a)5(i?'>o - a) 
which coincides with Eq. fll3ip up to the irrelevant constant det ~^(2a). 
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Our expression of the kinetic energy of the chain is shghtly more general, since both ends of 
the chain are free to move. 

An exhaustive discussion about the passage from discrete to continuous random chain systems 
may be found in [IQ] . 

Here and in the fohowing the interactions among particles of mass m mediated by the motion 
of the fluid are called hydrodynamic interactions, according to the definition of hydrodynamic 
interactions given in pj], Section 3.6, page 66. 



